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ABSTRACT 

Since pulsating subdwarf B (sdBV or EC14026) stars were first discovered (Kilkenny 
et al, 1997), observational efforts have tried to realize their potential for constraining 
the interior physics of extreme horizontal branch (EHB) stars. Difficulties encountered 
along the way include uncertain mode identifications and a lack of stable pulsation 
mode properties. Here we report on Feige 48, an sdBV star for which follow-up obser- 
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vations have been obtained spanning more than four years, which shows some stable 
pulsation modes. 

Wc resolve the temporal spectrum into five stable pulsation periods in the range 
340 to 380 seconds with amplitudes less than 1%, and two additional periods that 
appear in one dataset each. The three largest amplitude periodicities are nearly equally 
spaced, and wc explore the consequences of identifying them as a rotationally split 
i = 1 triplet by consulting with a representative stellar model. 

The general stability of the pulsation amplitudes and phases allows us to use the 
pulsation phases to constrain the timescalc of evolution for this sdBV star. Addition- 
ally, we are able to place interesting limits on any stellar or planetary companion to 
Feige 48. 

Key words: 

Stars: oscillations - stars: variables - stars: individual (Feige 48) 



1 INTRODUCTION 

To date, over 30 pulsating subdwarf B (EC 14026 or sdBV) stars have boon identified, with pulsation periods ranging from 68 
to 528 seconds and with amplitudes generally less than 50 millimagnitudes (mmag). For recent reviews of this class of stars, 
see Kilkenny (2001), and Reed, Kawaler & Kleinman (2000), for observational properties; Charpinet, Fontaine & Brassard 
(2001 and references therein) describe in detail some important aspects of pulsation theory in sdB stars. Most sdBV stars 
show periods at the short end of the range, and probably represent stars close to the Zero Age Horizontal Branch (ZAHB). 
PG 1605+072 is the longest period, and lowest gravity, sdBV star, with Feige 48 being an intermediate object. In general, 
the longer-period sdBV stars represent more highly evolved objects. 

Feige 48 was identified as a "faint blue star" as part of the Feige survey (Feige, 1958). It was re-categorized as an sdB star 
when it was observed as part of the Palomar-Green survey (Green, Schmidt, & Liebert, 1984). Koen et al. (1998; hereafter 
K98) identified five pulsation periods in Feige 48 in six observing runs from 1997 May to 1998 February. The periods detected 
by K98 range from 342 to 379 seconds with the largest amplitude being 6.4 mmag. Amplitude variability led K98 to conclude 
that mode beating was probably present implying that other unresolved modes were present in their data. This provided the 
motivation for our follow-up observations. Heber, Reid, & Werner (2000; hereafter HRW) obtained a high resolution (0.09A) 
spectrum of Feige 48, from which they determined log (?=5.50±0.05 and Teff=29500±300 K. This places Feige 48 among the 
coolest sdBV stars known with a surface gravity intermediate between PG 1605+072 and the rest of the class. 

Here we report on our multi-year campaign of high-speed photometry of Feige 48. In Section 2, we outline our observations. 
Section 3 describes the time series analysis and period identifications. We report on a stellar model fit to Feige 48 in Section 4. 
The phase stability of pulsations is described in Section 5, where we use this stability to place interesting limits on any possible 
planetary companion. Section 6 gives our conclusions and outlines future observations for Feige 48. 



2 HIGH SPEED PHOTOMETRY 

We began observing Feige 48 in 1998 November, with our most recent observations acquired in 2002 May. Table 1 provides a 
complete list of our observations (including K98's observations). We acquired most of the data using 3 channel photoelectric 
photometers as described in Kleinman, Nather, & Phillips (1996). At Fick Observatory, we used a 2 channel photometer 
of similar design. The Fick data typically have an ^30 minute gap during the night, as the telescope mount requires the 
photometer to be disconnected when the telescope exchanges sides on the pier. Because this instrument is a 2 channel 
photometer, the data were occasionally interrupted to measure the sky background. All photometers used Hamamatsu R647 
photomultiplier tubes. Data acquired at Calar Alto and SARA were obtained using CCDs with 5 second exposures on a 30 and 
15 second duty cycle, respectively. As both the target and comparison star were in the same CCD field, differential photometry 
removed extinction and sky variation. However, since extinction is wavelength-dependent, colour differences between the stars 
produced small non-linear trends in the data. We remove these trends by dividing by a low (2 - 4) order polynomial fit to the 
single-night data. Bad points in the CCD data were removed by hand. No filters were used during any of the photoelectric 
observations to maximize the photon count rate, whereas the CCD observations used filters to approximate the passband of 
the Hamamatsu phototubes (Kanaan et al. 2000). 

As Table 1 indicates, we have obtained a total of «380 hours of time series photometry on Feige 48. The data span from 
1998 January to 2002 May (the two runs in 1997 were too short and too temporally separated to be useful for this analysis). 
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Table 1. Observations of Feige 48 



Run Length Date Observatory Run Length Date Observatory 
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3 THE PULSATION SPECTRUM OF FEIGE 48 

We follow the standard procedure for determining pulsation frequencies from time series photometry obtained using the Whole 
Earth Telescope (see, for example, O'Brien et al. 1998). In short, we identify the principal pulsation frequencies with a Fourier 
transform of light curves of individual nights. We then combine data from several contiguous nights to refine the frequencies. 
Once the main frequencies are found, we then do a nonlinear least squares fit for the frequencies, amplitudes, and phases of 
all identified peaks, along with their uncertainties. 

To work around monthly and annual gaps between observing runs, we begin our analysis of the data in separate, relatively 
contiguous subgroups. The dates and data hours obtained for the subgroups are given in Table 2, the temporal spectra and 
window functions for the groups are plotted in Fig. 1. All groups were analysed independently, without using periods detected 
in other groups. This independent group reductions decrease the likelihood of selecting a daily alias over a real pulsation. 
Only in our Group X data, do we detect a mode (/2) inconsistent with the other group reductions. As such, we presume that 
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Table 2. Subgroups used in pulsation analysis. 



Group Inclusive dates Hours of data. 



I 


1998.28.01 - 


1998.01.02 


26.7 


II 


1998.22.11 - 


1998.26.11 


10.2 


III 


1999.06.03 - 


1999.13.04 


63.6 


IV 


1999.10.12 - 


1999.16.12 


17.5 


V 


2000.08.02 - 


2000.05.03 


42.4 


VI 


2000.04.05 - 


2000.05.05 


9.8 


VII 


2000.02.11 - 


2000.22.12 


45.7 


VIII 


2001.18.01 - 


2001.01.02 


55.0 


IX 


2001.19.04 - 


2001.30.04 


52.0 


X 


2002.05.04 - 


2002.21.05 


56.8 



Table 3. Comparison of the pulsation frequencies (in //Hz) detected in various runs. Formal least squares errors are provided in 
parentheses. 



Group 


/I 


/2 


/3 


/4 


/5 


/iHz 


It 


2636.96(15) 




2850.530(40) 


2874.40(23) 2877.310(50) 


2917.700(70)* 


III 


2641.98(1) 


2837.53(1) 


2850.833(4) 


2877.157(3) 


2906.275(4) 


V 


2641.49(2) 


2837.53(1) 


2850.833(3) 


2877.177(4) 


2890.025(19) 2906.299(8) 


VIII 


2642.00(7) 


2837.53(5) 


2850.818(17) 


2877.153(12) 


2906.266(22) 


IX 


2641.98(6) 


2837.78(3) 


2841.151(9) 2850.916(26) 


2877.185(13) 


2906.141(23) 


X 


2()il.S(ii(i) 


282(),!)7(())- 


28r)0. 811(12) 


287V.22(){1()) 


2!)(K).()()5( 13) 



t These frequencies are directly from K98. 

* Indicates modes offset by approximately the daily alias (11.56 fjMz). 



our periods for /I — /5 are not aliases, with the exception of /2 in Group X, which is a daily alias away from the real period. 
Frequencies determined for the bettor data sets arc in Table 3 with the corresponding temporal spectra of the best groups 
(Groups III, V, and IX) and the prewhitened residuals in Fig. 2. Though some signal remains in the Fourier transform after 
prewhitening within these groups, we cannot distinguish any remaining pulsation frequencies from noise. 

From our best solution for the Group III data, temporally adjacent groups were added one at a time and froquoncios 
and amplitudes were fit again until a satisfactory fit weis determined for all the data. Though there is still a chance that 
some modes may be off by an annual or monthly alias, the periods, frequencies and amplitudes in Table 4 represent our best 
solution. Also indicated by Table 3 are throe pulsation periods that appear only within the span of a single group. However, 
the amplitudes are sufficiently low that some possibility exists that they could be due to noise and or aliasing. As such, we will 
not include them in our analysis, but we mention them as they are possibly stochastically excited modes within an otherwise 
stable pulsation spectrum. Note also that an error corresponding to an annual (or even monthly) alias produces only a small 
change in the period or frequency, and so will not affect the modeling for asteroseismic analysis. Such a mistake, however, 
would be fatal to the period stability analysis. 

Fcige 48 shows five distinct and consistent pulsation modes. Four of the stable modes cluster with periods near 350 
seconds and a single mode lies at 378.5 seconds. The three shortest period modes also have the highest amplitudes - over 
three times higher than the two longer period modes. 

Our initial interest in follow-up observations of Feigc 48 was amplitude variability observed by K98. Our hope was to 
detect unresolved pulsations that could be responsible for the apparent amplitude variability they reported. However, it 
appears that the pulsations are resolved, and each has a variable amplitude of at least 30 %. The amplitudes determined for 

Table 4. Least squares solution to the entire data set. Formal least squares errors are in parentheses. 



Mode 


Period 


Frequency 


Amplitude 




(sec) 


(mHz) 


(mma) 


/I 


378.502960(37) 


2641.98731(23) 


1.19(6) 


/2 


352.409515(32) 


2837.60791(21) 


1..30(6) 


/3 


350.758148(3) 


2850.96729(7) 


4.17(6) 


/4 


347.565033(19) 


2877.15942(4) 


6.35(6) 


/5 


344.082794(15) 


2906.27734(8) 


3.57(6) 
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Figure 1. Temporal spectra and Window functions of Feige 48 for the groups listed in Table 2. 



each of our groups for the three high amplitude modes are shown in Figure 3. All 3 amplitudes show change, though only /3 
has a dramatic degree of variability. 



4 ANALYSIS OF THE PULSATION SPECTRUM OF FEIGE 48 
4.1 Frequency splittings in Feige 48 

For most sdBV stars, sufficient data do not exist to resolve the complete pulsation structure. For stars with resolved temporal 
spectra (Kilkenny 2001) there are typically too many modes to be accounted for by current pulsation theory unless high i 
values are included (where, by high £, we mean £ > 3). Though higher £ modes could be present, such modes sufTer from 
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Figure 2. Temporal spectra (top) and residuals after prewhitening by frequencies in Table 3 (bottom). 

severe cancellation effects across the unresolved stellar disk. In any case, identification of rotationally split multiplets {i = 1 

triplets or £ = 2 quintuplets of nearly equally spaced modes, for example) could aid in accounting for the many modes 
seen. Unfortunately, previous studies have been limited by a lack of equally spaced (in frequency or period) pulsations as a 
constraint on the observed £ values (with the exception of PG 1605+072; Kawaler 1999). 

Even though Feigc 48 has a relatively simple pulsation spectrum, understanding why this star pulsates with these 
frequencies still presents a problem. The tight cluster of four modes with periods spanning a range of less than 10 seconds is 
impossible to accommodate with purely radial pulsation modes in sdB models. Even appealing to nonradial pulsations, the 
closeness of these periods means that rotation (or other departures from spherical symmetry) must play a role. The reason 
is that if all are m=0 modes, standard evolutionary models of sdB stars at the same Tea and log y as Feige 48 do not have 
dense enough frequency spectra to account for these 4 modes, even if one appeals to modes of higher degree than £ = 3. 

However, Feige 48 may provide important clues in its observed pulsation spectrum. The frequency difference between /3 
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Figure 3. Pulsation amplitudes for the 3 largest amplitude modes. 



and /4 (26.2 fjHz) is very close to the diflerence between /4 and /5 (29.1 /iHz). Though not exact, this sphtting is highly 
suggestive that /3, /4, and /5 are a rotationally split mode, and probably £ = 1. We note that the frequency splitting is not 
precisely equal (26 versus 29 /uHz) . First order pulsation theory (if applicable here) says that the splitting should be precisely 
equal - though observations of rotational splitting in white dwarf stars show asymmetries such as this (i.e. in PG 2131+066, 
Kawaler et al. 1995 or GD358, Winget et al. 1994). Such asymetries can be caused by many higher-order processes including 
rotation and magnetic fields However this mainly observational paper is not the medium for such a discussion, so we will 
ignore the small departures from symmetry. 

Another possible clue is that the spacing between /2 and /3 is 13.2 /uHz, which is almost exactly half the frequency 
spacing between /3 and /4. Thus we could interpret the observed frequencies as follows: the /3, /4, /5 set make an £ = 1 
triplet, meaning the model needs to fit /4 with an £ — 1, and /3 and /5 reveal the rotation rate. Or, /2, /3, and /4 arc three 
components of an ^ = 2 quintuplet with a spacing of 13.2 /xHz. In this case the m = component could either be /3 itself or 
/3 + 13.2 /iHz (with a period of 349.1 seconds). 



4.2 Comparison with standard evolutionary models 

Though we by no means have a complete grid of models and cannot quantify the uniqueness of our results, we can see if 
either of the above possibilities is consistent with expectations from standard evolutionary models of sdB stars. We take the 
approach of Kawaler (1999); create several series of evolutionary models that pass through the spectroscopic errorbox of T^b 
and log 5, being sure to sample several regions. Prom these grids, we search for periods near those observed, creating a list 
of models and pulsation periods. Prom within this list, we further constrain the match using the required value of £ and 
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m for each possible interpretation of tlie splittings noted above. For the best model series, we create models with smaller 
evolutionary steps, finding the one with the best fit to the periods. 

With 5 froquoncios available, there arc multiple ways to match observations with model frequencies depending on the 
assumed values of £, and m for each mode. The most obvious assumption to make is that modes /3, /4, and /5 form a 
rotationally split triplet with £=1, m=-l,0,+l. This choice has no "missing members" of the multiplet. With this assumption, 
the model that fits the spectroscopic data must show an £ = 1 mode at /4, and modes at fl and /2. Other choices for this 
triplet which do require that some components of the multiplet are unobserved are £=2, with Am = 2 (m=-2,0,+2) or Am = 1 
(e.g. m=-2,-l,0 or m=-l,0,+l). Similarly, the interpretation of /2, /3, and /4 being components of an £ = 2 multiplet, as 
described above, is viable. 

Without choosing a priori one of these interpretations, we examined evolutionary models within our preliminary grid 
as described below. We generated full evolutionary stellar models using a version of the ISUEVO stellar evolution program 
(Dehner 1996; Dehner & Kawaler 1995; Kawalcr 1999) that incorporates semiconvection in the core helium-burning phase. 
Our initial grid of evolutionary tracks and models spans the observed range of Tefi and logg for sdB stars with a core 
mass of O.47M0 and solar metallicity. We computed evolutionary tracks for models with hydrogen envelope masses ranging 
from to 0.00550 Mq. From this initial grid, we focused on model series whose evolutionary tracks pass within la of the 
spectroscopically determined values of Teu and logp (HRW). For these models, we then calculated their pulsation periods to 
see which (if any) had radial or nonradial pulsation periods near those observed, for any possible choice of £ and identification 
of the m = component of a rotationally split multiplet. For the sequence that most closely matched the observed periods, 
we iteratively produced models with smaller differences in If shell masses and smaller evolutionary time-steps. 

We did find a model that matched the spectroscopic constraints, and, with an ad hoc rotational splitting of 27.73/iHz 
(the average of the observed splittings) , could explain four of the five observed frequencies. The closest model fit came from a 
model in the evolutionary track shown in Fig. 4. This model does an excellent job in explaining /2 through /5 and requires the 
identification of /3, /4, and /5 as a rotationally split £ = 1 triplet. The model is fairly evolved (as we suspected) and has nearly 
exhausted its core helium, with only 0.74% (by mass) of the core composed of helium. Table 5 shows a comparison between 
the observed periods and best-fitting model periods, and the observed spectroscopic properties and the model parameters. 
This model fits all but the lowest frequency to a precision of better than 0.2 s. The model temperature agrees to within 150 K, 
and log g agrees to within 0.02 dex of the measured values; these are well within the la uncertainties quoted by HRW. 

Despite the perception that there are many degrees of freedom, this procedure produced only 1 model series that fit the 
observations: i.e. no models had appropriate £=2 or £=1, n=2 periods. Of course, as we do not have model grids covering all 
parameter space, it is possible another model, with perhaps a different core mass or metalicity may fit better or with altogether 
different £ identifications. This procedure also failed to produce a model that could explain all 5 observed frequencies in terms 
of normal modes and rotational splitting. Though the model does produce an £=3, m=0 mode at 374 s period, which is near 
the observed 378 s mode, wc currently find no evidence to suggest that Fcige 48 has £ > 2. As such, wc must confess that our 
model does not reproduce this period without appealing to high £. Without further evidence (such as other observed members 
of the multiplet) for invoking high £ modes, we are forced to leave fl as unmatched by our model. Additionally, any £=2 
matches fitted less observed periods. Since wc do not have a complete sample of models and this is really just an illustrative 
example, we are not alarmed. However, it could also indicate that our current models do not include enough physics to be 
accurate. 

The pulsation results for the closest-fitting modes in our best-fitting model series are shown in Fig. 4. Even a small change 
in log g (as an indication of age) of 0.005 dex changes the calculated periods by more than 3 seconds, worsening the fit to the 
observations. Likewise, a change in the envelope layer thickness quickly destroys the fit by moving the evolutionary track's 
path away from the spectroscopic error box. As the right panel of Fig. 4 shows, within this period space, model periods are 
relatively uncrowded. Overtones are separated by ~100 s for £=0 and 1 modes, and ~50 seconds for £=2 modes. Overtones 
for £=2 do appear in the top right and lower left of the figure. 

4.3 Testing the mode identifications 

A test of our (or any) model is the measured constraint on rotational velocity. The observed average rotational splitting 
of 27.7/iHz imposed on our £ = 1 model identification implies a rotation period of 0.42 days (10 hours). With a radius of 
0.20 R©, this model has an equatorial rotation velocity of 24 km s^^. To match the constraints of HRW (w sini < 5km ■ s~^) 
requires i < 12°. If we use the less restrictive value determined by HRW for only the unblended spectral lines in Feige 48 of 
vsini < 10km ■ s^^, our inclination limit increases to i < 25°. 

The identification of am £ = 1 triplet and a radial mode suggest an observation that can be used to test the model. As in 
white dwarfs (Kepler et al. 2000), time series spectroscopy (particularly in the ultraviolet) should present an effective means 
of discerning between low and high order {£) nonradial pulsations. Though it is still in its infancy for sdBV stars (O'Toole et 
al. 2002; Woolf et al. 2002), if the 378 s period is £=3, it should be obvious in UV spectroscopy (perhaps less so in the optical) 
where it should have a significantly higher amplitude than at optical wavelengths. The same is true for our identification of the 
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Table 5. Comparison of observations with our best-fitting evolutionary model. 
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Figure 4. Comparison of the model to the observations. The left panel shows the evolutionary model track (solid line) containing our 
best-fitting model (dot). The dashed line is the ZAHB and the rectangle is the spectroscopic Icr error box. The right panel compares the 
model periods (points) to those observed (solid lines). The vertical dashed line indicates the best-fitting model. 



£=1 triplet. If any member of our identified triplet really has a different £ value, the wavelength dependence of its amplitude 
will be different. Such a test should be obtained as an independent confirmation of our £ = 1 determination. While this test 
can be applied to any sdBV star, Feige 48 has comparatively long pulsation periods (exceeded only by PG 1605-1-072) and 
its rather simple temporal spectrum (only 5 periods compared to 55 for PG 1605-1-072) make it an ideal candidate for time 
series spectroscopic study. 



5 STABILITY OF THE PULSATION PERIODS 

As a star evolves, the pulsation properties evolve in response. In the case of the subdwarf B stars, evolutionary models indicate 
that they reside on or near the ZAHB for approximately 10* years. Upon exhausting their core helium supply, they leave the 
HB, their \ogg goes down, and pulsation periods lengthen. The time scale for evolutionary pulsation period changes (P/P) 
after leaving the HB is about 10 times faster than while on the HB. 
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If Feige 48 has a comparatively small log g because it is a mature HB star that has left the ZAHB, we expect Feige 48 to 
have an evolutionary P smaller than PG 1605+072, yet larger than for shorter period pulsators. Since PG 1605+072 does not 
appear to have pulsations stable enough for an analysis of secular period change caused by evolution (Reed 2001), Feige 48 
is the best candidate to examine the e-folding time for structural changes caused by its core evolution. As a guide, the model 
described in Section 4 has a P=1.714 • 10~® syr~^. With ~3 years of usable data, the phase of a 350s period should change 
by ~14 seconds in that time. This is close to our limit of detection. 

To examine long-term phase change, we followed methods outlined in Winget et al. (1985) and Costa & Kepler (2000). 
First, we obtain a best-fitting least squares fit to all of the periodicities present over the entire span of the observations (Table 
4). We then fix the frequencies at these best-fitting values, and recompute the pulsation phases (again via least squares) for 
each group in Table 2 (note that Groups III and V were divided into two subgroups each because of the long length of the 
runs). This computed phase represents the observed time of maximum (O) for that group, which differs from the computed 
time (C) from the fit to all data. The resulting O ~ C diagram is shown in Fig. 5 for the three highest amplitude modes (with 
the pulsation period indicated in each panel). The phase zero point is that defined in K98 as JD=2450571.50. 

Fig. 5 shows that the phases are stable throughout most of our observations. Until the Group X data were collected, 
we believed the Group I data suffered a timing error (proposed in K98). However, it now seems apparent the pulsation 
modes were only stable over a limited timespan (from Group II through IX). As a check, we analyzed various subgroups of 
Group X data and reproduced the same phase result. Such a problem is observed in other sdBV stars over a much shorter 
timescale (Reed 2001). Though we are disappointed in the apparent lack of phase stability, we still have an approximately 
three year span of phase-stable observations. We therefore used the phase-stable data to place upper limits on the magnitude 
of P/P. The combined O — C data of the three highest amplitude modes were weighted and fitted using least squares and 
give P/P=4.9 ± 5.3 x lO^^^s^^. The data are therefore consistent with zero period change and provide a la lower limit on 
an evolutionary timescale of 3.1 x 10*^ years. Of course evolution is not the only thing that can drive period changes (see for 
example Paparo et al. 1998). However, evolutionary models predict that sdBV periods should change in a predictable way 
(increasing just off the ZAHB, then decreasing after core helium exhaustion, and finally increasing again during shell helium 
fusion). By measuring P for several sdBV stars at different stages, we should be able to determine if evolution (as predicted) 
is driving the period changes. 

5.1 O — C variations from reflex orbital motion - planets around Feige 48? 

With the phase-stable portion of the data, we can also place useful limits on companions to Feige 48. Orbital reflex motion 
would create a periodic shift in the pulsation's arrival time, which would be observed in pulsation phase. Thus any companion 
must create a periodic phase change within our O — C uncertainties over a scale of days to years*. To place limits on 
companions, we calculated companion mass as a function of binary period and semimajor axis by fitting sine curves (for 
circular orbits) within our 1, 2, and 3(7 O — C limits. To assure an upper limit (Msini) on reflex motion, we use the "noise" in 
the FT of our O — C as a Icr lower limit: This results in a minimum phase shift of 5 seconds for binary periods under 20 days 
and 4 seconds for longer periods. Fig. 6 graphically presents our sensitivity to orbital companions. The top line represents 
the 3(7 limit for i = 12° (our model constraint). In the short period case (periods under 30 days), the constraint is the limit 
imposed by the phase errors of individual runs within the data; in the long period case it is the flatness of the (O — C) diagram 
(including the errors of combined runs) over the phase stable region of our observations. The drop at 30 days corresponds to 
the change from O — C values determined for single runs to group data sets. 

Feige 48 is a horizontal branch star that has lost considerable mass between the red giant branch and its current evo- 
lutionary state. Any companion separated by more than ~1 AU is far enough away that common-envelope evolution (or 
vaporisation) has been avoided. In addition, the orbital separation will have roughly doubled as Feige 48 lost approximately 
half its mass during the red giant phase. This should produce two cases for binaries: Close stellar companions with original 
separations < lAU will produce a short period binary (which have periods on order of weeks or less) after a common envelope 
phase. Companions distant enough to avoid a common envelope phase (or vaporisation) will have orbital periods on the order 
of a year or more. The two panels of Fig. 6 reflect this duality (though it does cover all periods between 2 days and 5 years). 

The left panel of Fig. 6 indicates our limits on stellar companions. Our la limit is less than O.IM0 sini for a binary period 
of 3 days. The right panel shows our limits on sub-stellax companions. Our "average" 3a limit for i = 12° is w 12M^, while 
our best la limit would detect Jupiter at a period of 2.5 years. Our data are currently sensitive to the extra solar "warm 
Jupiter" type planets being detected^ at a distance of 0.6 - 3.0 AU. Planets with orbital separations less than ^1 AU would 
not have survived the red giant phase. Our data do not rule out a companion in an extremely short period binary or at low 
inclination. 

* We assume we could detect an orbital period up to twice our observed time base. 

t A complete list of extra solar planets is maintained at http://www.obspm.fr/encycl/catalog.html. 
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Figure 5. O-C diagrams for the 3 largest amplitude modes. 



6 CONCLUSIONS 

From our multi-season photometry of Feige 48, we have consistently detected five pulsation periods. Of these five, three (/I, 
/3 and /4) are consistent with K98. One frequency (/5) differs by a daily alias, while K98's fifth frequency (2874 fjHz) is 
not detected in our data. In data sets V and IX, we also detect new pulsation frequencies at 2890 and 2843 /iHz respectively, 
indicating that there may be some stochastically excited pulsations in Feige 48. This is consistent with pulsation behavior 
seen in another sdBV star, PG 1605+072 (Reed 2001). 

Our attempted model fit follows the strategy that has been successfully applied to other classes of pulsating stars, but has 
rarely worked for sdBV stars; namely using observed frequency splittings to impose £ constraints on models. Using standard 
evolutionary models, our preliminary model grid includes a model that is able to explain all but the lowest frequency stable 
pulsations. Most exciting is the fact that such standard models fail to explain all five frequencies. Thus, despite its relative 
simplicity and the richness of the parameters available, the failure of this model suggests that standard stellar evolution theory 
does not fully explain the evolution of sdB stars and or the nature of pulsations within them. We have something new to 
learn. 
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Figure 6. 1<t, 2a (dashed line), and 3it upper limits on companions to Feige 48. The time axis is continuous between panels, but changes 
scale from days to years. The mass axis is discontinuous between panels as the left panel is in solar masses and the right panel has units 
of Jupiter masses. 



Our modeling example shows that Feige 48 should also serve as an interesting test for other methods of mode identification. 
Though optical multicolour photometry was not useful for identifying pulsation modes in KPD2109+4401(Koen 1998), we 
expect that UV multicolour photometry as developed by Robinson, ct al. (1995), will be useful to determine i/ high-order £ 
modes are present in sdBV stars (eis indicated by Brassard et al. 2001 and Billeres et al. 2000). The argument for high-order 
I values is particularly interesting in light of the frequencies detected in Group V's data. If the lowest frequency mode is 
disregarded, the remaining modes have frequency spacings of 13.3, 26.4, 12.8, and 16.3 /iHz respectively. If these were all parts 
of a single, rotationally split mode, it would require £ > 3. Such an I value should be apparent in UV multicolour photometry 
(see, for example, Kepler et al 2000). As such, we look forward to the analysis of HST data obtained by Heber (2002). Should 
Heber's (2002) HST data agree with our £=1 interpretation, Feige 48 would make an excellent star to calibrate other mode 
identification methods in sdBV stars such as optical time-series spectroscopy (O'Toole et al 2002; Woolf et al 2002). 

The results of our O — C analysis are consistent with a non-binary nature for the star within the data limits. It also 
indicates that using the O — C diagram to detect planets around evolved stars is possible, though in this case we did not detect 
any. We plan to continue to monitor Feige 48 over the next several years to tighten the constraints on planetary companions. 

Our limit on a stellar companion also addresses the origin of sdB stars. Binary evolution is a candidate for producing 
sdB stars, either through common-envelope evolution (Sandquist, Taam & Burkert 2000; Green, Liebert & Saffer 2000) or 
via Roche-lobe overfiow near the tip of the red giant branch (Green, Liebert & Saffer 2000). Though observations indicate 
that a great many sdB stars are in binaries (Reed & Steining 2003; Green, Liebert & Saffer 2000; Han et al 2002), either 
the evolutionary sequence that produces sdB stars is independent of binary evolution (D'Cruz et al 1996), is bimodal, or has 
several paths that can result in the production of an sdB star (perhaps including the merger of two low mass white dwarfs as 
described by Iben & Tutukov 1986). For the case of Feige 48, it would appear that it is either in a short period binary (whose 
orbital period is commensurate with the ^^10 hour rotation period predicted with our model), in a long period binary with an 
orbital period substantially longer than our data (which would rule out Roche-lobe overfiow, so the companion would have 
no effect on the evolution of the pre-sdB star), in a binary with an extremely low inclination or a single star. 
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